Non-Markovian hole spin kinetics in p-type GaAs quantum wells 
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Based on fully microscopic kinetic spin Bloch equation approach, we show the non-Markovian 
effect in spin dephasing of heavy holes in p-type GaAs quantum wells. The non-Markovian effect 
manifests itself in spin dephasing when the mean spin precession time is shorter than the momentum 
relaxation time. The spin dephasing becomes slower when the non-Markovian effect is important. 
Moreover, quantum spin beats due to the memory effect of the non-Markovian hole-longitudinal 
optical phonon scattering are predicted. 

PACS numbers: 72.25.Rb, 78.90.+t, 71.10.-w 



Much attention has been devoted to the spin de- 
gree of freedom of carriers in semiconductors re- 
cently for the purpose of the possible application 
of spintronic devices ,L2i£i4 Understanding spin relax- 
ation/dephasing is one of the most important issues^ 
Recently an extensive investigation has been performed 
based on kinetic spin Bloch equations (KSBEs)^ to un- 
derstand the spin relaxation/dephasing in zinc-blende 
semiconductors /' 8 ' 9 ' 10 ' 1, ' In these works, together with 
earlier studies^ the kinetics was treated in Markovian 
limit. The Markovian limit is a good approximation 
in the strong scattering limit where the spin precession 
time is long compared to the momentum relax- 

ation time t p , i.e., fl^ 1 3> t p , which is mostly the case 
for n-type GaAs quantum wells (QWs)2 unless at very 
low temperature. 12 However, for hole systems (especially 
heavy hole systems) the situation changes a lot due to 
the strong spin-orbit coupling strength. Here the mean 
spin precession time becomes comparable with the mo- 
mentum relaxation time, i.e., fl^ 1 < t p . Then Marko- 
vian approximation is not appropriate any more. In this 
work we are going to extend our previous kinetic spin 
Bloch equation approach^ to the non-Markovian limit to 
investigate the spin dephasing of heavy holes in p-type 
GaAs QWs. In non-Markovian kinetics, the energy con- 
servation of the scattering is lifted and there is memory 
effect due to the history dependence of scattering. In 
fact, the non-Markovian effect has been extensively ex- 
plored in the ultrashort-pulse spectroscopy in semicon- 
ductor optics / 3 : 14 Recently it has been explored in quan- 
tum dot system for spin dephasing induced by hyperfine 
interaction*^ Glazov and Sherman also investigated the 
electron spin relaxation in QWs with strong magnetic 
field by Monte Carlo simulation^ 

We start our investigation from a p-type GaAs (100) 
QW with well width a being small enough so that 
the heavy-hole and the light-hole bands can be treated 
separately i 1 1 1 1 7 i 18 We only consider the heavy hole in this 
report. Moreover, only the lowest subband is considered. 
The KSBEs constructed by the nonequilibrium Green 



function method read^* 9 - 



Pk — Pk\coh + Pk\scat 



(1) 



in which pk represent the density matrices of heavy hole 
with momentum k. The diagonal terms pk,<ra = fka 
(a = ±3/2) represent the hole distribution functions 
and the off-diagonal ones p t 3 3 = p* 3 3 describe the 

rK *>2 2 k ,~2 2 

inter-spin-band correlations for the spin coherence. The 
coherent terms pk\coh describe the coherent spin preces- 
sions around the effective magnetic field from the Rashba 
terms fl(k) and its expressions can be found in Ref. [Til . 
In the present report, we only consider the heavy hole- 
longitudinal optical (LO) phonon scattering in pk\scat 
which is written in the non-Markovian limit 



dpk 

dt 
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sum , 
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with 



dr[N> P <{t)p>{t) 
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Here p£ = p k and p^ 
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Pk- 

T 



N and N' 



N< 

N a + 1 with N = (e^o/fcB'J _ ly 1 standing for the 
Bose distribution of the LO phonons with frequency u>q 
at temperature T. = h 2 k 2 /2m* is the energy of hole 
with wave vector k and effective mass m* . gk-k',?, is the 
hole- phonon interaction matrix element.— 

It is noted that the KSBEs in non-Markovian limit 
are integro-differential equations. The integral terms in 
indicate the history dependence of scattering 

scat 

processes which lead to the memory effect, with respect 
to that all the density matrices enter only at the earlier 
time r. However, if p^ change slowly on the time scale of 
collision, one can disregard the retardation and pull the 
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distributions out of the time integral. In this completed- 
collision approximation, one can simply get in the scat- 
tering term ^-functions denoting energy-conservation. In 
this way one comes to the Markovian limit, with memory 
effect lifted. For comparison, the KSBEs will be solved 
in both Markovian and non-Markovian limits. 

The integro-differential KSBEs in non-Markovian limit 
can be transformed into a larger set of differential equa- 
tions by factorizing the integral terms and making time 
derivative actions on them.— Finally the group of equa- 
tions to be solved numerically are written into 



dt 



dt 



k'<j 2 



(4) 



kk' 



dt 



i[uj Q - (E k > - Bk)/ft]5kk' 

+ [N > P^(t) P >(t)~N<p>(t)p<(t)} . (5) 



Once the KSBEs are solved numerically, the temporal 
evolution of the hole distribution /k±s(t) and the spin 
coherence yO k 3_|(t) are obtained. As discussed in the 
the spin dcphasing can be ob- 



,19,20 



previous papers j* 
tained from the slope of the envelop of the incoherently 
summed spin coherence puH{t) = Sk |Pk|-|(*)l- 

We solve the KSBEs numerically and the main re- 
sults are plotted in Figs. [1] and O In the calculation, 
T = 300 K and a = 5 nm. The hole density N = 
N + + A_ = 4 x 10 15 na- 
tion P = [N + -N_)/(N+ 
standing for the density of holes with up- (down-) spin. 
The material parameters of GaAs QWs can be found in 
Ref. [ll| except that the heavy hole effective mass m* , de- 
duced from the 8x8 Luttinger-Kohn Hamiltonian based 
on the quasi-degenerate theoryj2i depends on the width 
of the QW and is 0.237mo when a — 5 nm. We select 
such a value of E z that it satisfies r )l ] l 7h E z m /h 2 = 1.25 
nm in the Rashba term when a = 5 nm*ii 



■ 2 and the initial spin polariza- 
+ N_) = 0.025 with N + (iV_) 
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FIG. 1: (Color online) Time evolution of the incoherently 
summed spin coherence phh- Solid curve: in non-Markovian 
limit; Dashed curve: in Markovian limit. 



In Fig. [T] the time evolutions of the incoherently 
summed spin coherence in both the Markovian and non- 
Markovian limits are plotted. One can see that the spin 
dephasing is slower in the non-Markovian limit than the 
Markovian one. This is in agreement with the result of 
Glazov and Sherman^ and can be interpreted as the in- 
fluence of the memory effect — in non-Markovian kinet- 
ics, the spin coherence of heavy holes can be partially 
kept due to uncompleted hole-LO phonon scattering pro- 
cesses. It is interesting to note that there are quan- 
tum spin beats superimposed in the decaying signal in 
the non-Markovian spin kinetics with the beat period 
being the period of the LO phonon (2itu)q w 0.115 
ps). This fascinating phenomenon is deemed as the char- 
acteristic feature of the non-Markovian effect, similar 
to the quantum beats in the ultrafast four-wave mix- 
ing signals in GaAs . 13 ' 14 The beating is caused by the 
transfer of spin coherence of holes between uncompleted 
scattering — when a phonon is partially emitted, it can 
be absorbed back by the hole system without losing co- 
herence. One can realize that this can not happen in the 
Markovian kinetics, where all the scattering processes are 
instantaneous and complete without feeding back. 
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FIG. 2: (Color online) Time evolutions of the incoherently 
summed spin coherence phh at scaling parameters \ — 5 (a); 
2 (b); 0.5 (c) and 0.2 (d), respectively. Correspondingly, the 
mean spin precession time fi _1 is 0.046 ps (a); 0.116 ps (b); 
0.463 ps (c) and 1.160 ps (d). Solid curve: in non-Markovian 
limit; Dashed curve: in Markovian limit. 



In order to gain more insight into the condition for non- 
Markovian effect, we repeat the calculations with differ- 
ent mean spin precession times. To do this, we introduce 
a scaling parameter x m t ne Rashba term, i.e., x^k- 
Experimentally the value of the Rashba coefficient can 
be tuned by bias voltagei 22 ' 23 The results are plotted in 
Fig. [U It is shown from the figure that the smaller the 
mean precession time f2 _1 = 2irh(xQ] s )~ 1 is, the more 
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pronounced the non-Markovian effect becomes. Hen 
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When x = 0.2 for which f2 _1 is much longer than the mo- 
mentum relaxation time T p) 2i the non-Markovian curve 
approaches the Markovian one (f2 _1 is about 0.232 ps 
when % = 1 and r p is about 0.121 ps). This suggests 
that when 3> r p , the Markovian approximation is a 
good approximation, as expected. However, when f2 _1 
decreases, the non-Markovian curves deviate from the 
Markovian ones and become more pronounced for smaller 
f2 _1 , indicating stronger non-Markovian effect. More- 
over, it is noted from the figure that strong quantum 
spin beats appear when the mean spin precession time is 
comparable with the momentum relaxation time. When 
the mean spin precession time is too short, i.e., the inho- 
mogeneous broadening 6 is too strong, the quantum spin 
beats are smeared out as shown in Fig. [2] (a) when x = 5. 



In conclusion, we investigate the non-Markovian spin 
kinetics of heavy holes in p-type GaAs QWs. The result 
is compared with the Markovian spin kinetics. We show 
that the non-Markovian effect — slower spin dephasing 
appears when the inhomogeneous broadening is strong 
(i.e., the mean spin precession time is shorter than the 
momentum relaxation time). We further predict quan- 
tum spin beats due to the memory effect of the non- 
Markovian hole-LO phonon scattering in the spin signals. 
Strong quantum beats appear when the mean spin pre- 
cession time is comparable to the momentum relaxation 
time. 
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